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LONG-TERM  GOALS 

The  long  term  goals  of  this  collaborative  research  effort,  which  is  coordinated  with  other  ONR  and 
Navy-funded  biological  research  projects,  are  two-fold.  We  aim  to  first  determine  how  the  distribution, 
abundance,  and  behavior  of  prey  affects  the  foraging  behavior  and  ecology  of  baleen  whales  off  the 
California  coast.  Baleen  whales  employ  a  variety  of  feeding  strategies  that  relate  to  the  behavior  of 
their  prey  and  understanding  these  is  paramount  to  being  able  to  assess  changes  in  their  feeding 
ecology  arising  from  a  host  of  natural  and  human  factors.  Second,  these  empirical  findings  are  directly 
applicable  and  are  extremely  useful  in  determining  how  prey  affects  whale  behavior.  Given  the  broader 
goals  of  investigating  potential  behavioral  responses  of  whales  to  mid-frequency  military  sonar,  the 
ability  to  more  completely  describe  and  quantify  behavioral  responses  of  baleen  whales  to  controlled 
exposure  experiments  while  including  the  effects  of  prey  provides  a  novel  and  powerful  insight  into 
interpreting  responses  to  sound  and  controlling  for  environmental  factors.  The  basic  measurements  of 
foraging  ecology  and  behavior  also  provides  a  critical  means  of  interpreting  potential  responses  buy 
describing  the  energetic  consequences  of  any  observed  changes  in  behavior.  In  order  to  determine 
whether  and  how  behavioral  changes  occurring  in  baleen  whales  during  controlled  exposure 
experiments  are  related  to  sound  in  their  environment,  we  need  to  better  understand  and  quantify 
whether  and  how  changes  in  their  prey  environment  account  for  the  behavioral  change  as  well.  In 
baleen  whales,  the  behavioral  states  most  commonly  observed  are  feeding,  traveling,  resting,  and 
socializing.  Blue  whales  visit  the  southern  California  Bight  in  the  summer  months  primarily  to  forage, 
and  therefore  understanding  baseline  behavior  (such  as  how  changes  in  their  prey  affect  the  likelihood 
of  changing  behavioral  states)  is  necessary  to  adequately  describe,  understand,  and  effectively  mitigate 
the  affects  of  anthropogenic  sound,  including  military  sonar,  on  these  animals. 
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OBJECTIVES 


The  overall  objectives  of  this  multi-year  study  are  to  obtain  empirical  synoptic  measurements  of  fine- 
scale  prey  distribution  and  whale  diving  and  foraging  behavior  in  order  to  better  understand  baleen 
whale  foraging  ecology  and  better  interpret  responses  to  experimental  sound  exposure.  The  current 
project  has  already  enabled  us  to  obtain  basic  distribution  and  density  information  for  prey  concurrent 
with  foraging  mysticete  cetaceans  during  tagging  with  fine-scale  movement  sensors  in  the  context  of 
behavioral  response  studies  (specifically  the  Southern  California  Behavioral  Response  Study,  or 
SOCAL-BRS).  The  new  results  presented  here  clearly  and  quite  powerfully  demonstrate  that  data  on 
the  distribution  and  abundance  of  prey  are  essential  in  fully  understanding  how  changes  in  whale 
behavior  related  to  the  presence  of  human  sounds  are  mediated  by  these  factors  and  associated 
environmental  variables.  Subsequent  CEEs  involving  potential  behavioral  changes  in  foraging  marine 
mammals  should  build  on  these  novel  techniques  and  measurements  in  order  to  fully  describe  potential 
responses  (or  lack  thereof)  to  sound  exposure. 

APPROACH 

SIMRAD  EK60  echosounder  units  (38  and  120  kHz  echosounders  and  GPTs)  and  top-side  hardware 
were  made  available  for  the  project  through  collaborations  with  research  partners  at  Duke  University 
(Dr.  Doug  Nowacek).  A  specialized  echosounder  mount  and  towfish  appropriate  for  the  SOCAL-BRS 
platform  was  fabricated  with  support  from  this  award.  The  two  Echosounder  units  are  shown  mounted 
on  the  towfish  in  the  figure  at  the  right.  The  smaller  orange  echosounder  is  the  120  kHz  unit  and  the 
larger  echosounder  is  the  38  kHz  unit. 


Figure  1.  Specialized  echosounder  mount  and  towfish  used  in  the  SOCAL-BRS. 

Several  computers  were  required,  including  a  ruggedized  laptop  computer  specific  to  data  acquisition 
(Dell  Latitude  E6420  ATG)  and  a  laptop  computer  for  field/lab  data  analysis  (MacBook  Pro  - 
15.4/2. 3/2X2 GB/750/SD/HR- AG).  Additionally,  several  external  data  storage  drives  and  other 
ancillary  gear  (e.g.,  handheld  GPS  unit)  were  required  and  obtained  with  this  award.  The  information 
from  each  echosounder  is  processed  in  a  GPT  (general  purpose  transceiver)  that  also  acts  as  a  power 
supply.  The  data  are  then  streamed  through  Ethernet  cables  to  the  laptop  where  they  are  processed  in 
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customized  visualization  and  analysis  software  (Echo view).  The  data  are  stored  directly  on  the  laptop 
computer  and  then  backed  up  on  multiple  external  hard  drives  routinely.  A  hand-held  Garmin  GPS 
unit  is  connected  to  the  laptop  to  provide  a  time  and  location  stamp  for  the  echosounder  data  as  it  is 
acquired. 


Figure  2.  Left  -  Net  assembly  being  deployed  from  the  SOCAL-BRS  research  vessel. 

Right  -  example  of  prey 


Finally,  a  plankton  sampling  system  was  obtained  for  the  second  leg  of  the  project,  given  the 
interesting  measurements  made  initially  and  the  inability  to  obtain  plankton  species  identifications. 
This  consisted  of  a  triple-stitched  plankton  net  (100  cm  X  500  cm  X  1000  Microns  x  1 1  cm  diam; 
4.5”0  COD  end  aperature),  an  SS-ring  and  bridle  assembly,  a  complete  2-PC  PVC  COD  end  assembly 
(1000  microns),  and  a  mechanical  flow-meter.  This  net  assembly  is  shown  being  deployed  from  the 
back  of  the  SOCAL-BRS  research  vessel  in  the  photo  to  the  left  with  a  sample  of  krill  on  the  right.  In 
order  to  quantify  the  density  and  biomass  of  prey  measured  from  the  echosounders,  we  need  to 
generate  length-frequency  estimates  of  the  actual  targets  that  are  being  measured.  Thus,  incorporating 
the  net  into  our  sampling  protocols  will  allow  for  more  accurate  and  quantitative  estimates  of  prey  for 
our  analysis. 

Fine-scale  prey  density  and  distribution  and  individual  predator  behavior  was  measured  in  two  phases 
in  SOCAL-1 1  (late-July  to  mid- August  and  September  2011)  and  in  the  first  phase  of  SOCAL-12 
(July- Aug  2012)  using  the  existing  research  platform  (R/V  Truth).  By  analyzing  prey  and  predator  at 
fine  scales  (100s  of  meters),  we  can  begin  to  test  for  the  relationships  between  prey  distribution  and 
predator  behavior  and  understand  the  ecological  decisions  made  by  individual  whales  when  foraging, 
and  how  the  broader  oceanographic  environment  affects  blue  whales  in  southern  California. 

Prey  sampling  -  Prey  distribution  and  abundance  was  continuously  measured  using  38  and  120  kHz 
SIMRAD  EK60  echosounders  at  fine  scales  (<10  km).  Acoustic  data  collected  in  the  absence  of 
sighted  or  tagged  whales  were  treated  as  a  control  measure  of  ambient  prey  density.  Fine  scale 
sampling  methods  are  dependent  on  the  behavior  of  the  tagged  whale  so  an  iterative  approach  to 
sampling  prey  is  employed.  If  the  tagged  whale  is  traveling  (>lkm  per  hour  displacement),  a  zig-zag 
design  was  used  to  survey  prey  distributions  passed  over  by  the  whale  by  sampling  in  its  wake 
(~1 .4km  long  transects).  When  focal  whales  were  surface  feeding  (defined  as  observing  the  animal 
with  its  mouth  gaped  or  bubbles  located  where  the  animal  surfaced),  a  clover  leaf  sampling  design 
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allowed  the  measurement  of  prey  abundance  and  distribution,  with  the  center  of  the  sampling  box 
centered  around  the  whale  (See  figure  to  right).  The  sampling  design  around  non-feeding  and  non¬ 
traveling  (i.e.  resting)  individuals  is  identical,  with  a  cloverleaf  used  to  examine  the  prey  distribution  in 
the  absence  of  feeding.  When  measuring  prey  relative  to  surfacing  events,  transects  were  designed  to 
pass  within  500  meters  of  the  tagged  whale.  Correlations  between  whale  behavior,  prey  data,  and 
environmental  data,  will  only  be  considered  in  analysis  within  a  500m  radius  of  a  whale  surfacing. 

This  will  allow  us  to  quantify  the  distribution,  abundance  and  dimensions  of  prey  patches  in  close 
proximity  to  foraging  and  non-foraging  whales.  We  will  also  compare  the  two  frequencies  of  acoustic 
data  to  differentiate  krill  from  larger  fish  targets  as  krill  have  greater  backscatter  at  120kHz  than 
38kHz. 


Figure  3.  Clover  leaf  sampling  design  around  tagged  whale.  Each  leaf  is  lkm  from  the  center. 

Whale  data  -  Whale  behavior  (e.g.,  feeding/non-feeding)  is  inferred  from  the  tag  record  in  combination 
with  near  continuous  daytime  focal  surface  observations.  Tags  were  attached  from  ~6m  rigid-hulled 
inflatable  boats  (RHIB)  by  taggers  using  hand-held  poles  from  which  Woods  Hole  Oceanographic 
Institution  (WHOI)  Digital  Acoustic  tags  (DTAG)  were  deployed.  The  DTAG  is  a  small,  lightweight, 
pressure  tolerant  tag  capable  of  recording  data  for  up  to  ~20  hours  and  attached  to  the  whale  via 
suction  cups.  The  DTAG  measures  the  acceleration  in  the  animal’s  pitch,  roll,  and  heading,  as  well  as 
depth,  and  water  temperature  at  50  Hz.  The  tags  also  measure  sound  and  calibrations  have  been  made 
between  vertical  acceleration  and  flow  noise  to  determine  when  whales  lunge  underwater.  This  is 
determined  by  increased  acceleration  as  the  whale  approaches  a  prey  patch  and  dramatic  deceleration 
when  the  animal  opens  it  mouth  to  lunge  and  engulf  prey.  This  approach  has  been  published  and 
ground-truthed  for  several  species  of  baleen  whales,  including  blue  whales  and  thus  is  considered  the 
most  accurate  way  of  determining  feeding  events  in  baleen  whales  from  tag-derived  records.  Data 
from  the  pitch  record  also  allows  for  analysis  of  fluke  stroke  rates  and  relative  stroke  amplitudes  and 
combined  with  behavioral  observation  allows  the  identification  of  surface  feeding  bouts  and 
quantification  of  their  duration.  All  sensor  data  are  stored  in  flash  memory  on  the  tag  and  are 
downloaded  via  an  infrared  connection  to  a  computer  for  analysis.  The  tag  has  a  VHF  antenna  that 
transmits  when  at  the  surface,  allowing  us  to  follow  the  whale  when  it  is  either  out  of  visual  range  or 
during  nighttime.  Focal  follows  were  conducted  from  RHIBs  such  that  animal’s  position  was  recorded 
by  marking  a  GPS  position  at  the  location  (foot-print)  and  time  where  the  tagged  whale  made  a 
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terminal  dive.  Additionally,  we  augmented  this  method  by  collecting  high-resolution  range  and 
bearing  measurements  using  a  laser  range-finder  (Leica  Vector  IV),  to  georeference  the  surfacing 
locations  of  the  tagged  whale  more  frequently.  Similar  to  previous  studies  using  non-linear 
generalized  additive  models,  in  analysis  we  will  quantify  the  effects  of  remotely  sensed  environmental 
features  and  prey  abundance  on  the  distribution  and  abundance  of  whales  at  the  seascape  scale.  This 
approach  will  provide  estimates  of  (1)  prey  and  environment  in  the  functional  study  area  around  blue 
whales  and  (2)  the  functional  relationships  between  prey  density  and  school  size  and  predator 
aggregation  size. 

WORK  COMPLETED 

•  Prey  mapping  data  were  collected  concurrent  to  baleen  whale  tagging  and  playback 
experiments  on  during  each  of  the  SOCAL-12  and  SOCAL-13  field  efforts  while  towing  the 
prey  mapping  towfish  (8  total  weeks  of  field  effort). 

•  Submission  of  the  manuscript:  “Feeding  performance  by  sympatric  blue  and  fin  whales 
exploiting  a  common  prey  resource”,  by  Ari  S.  Friedlaender,  J.  Goldbogen,  E.  Hazen,  J. 
Calambokidis,  and  B.  Southall,  to  Marine  Mammal  Science. 

•  Analysis  of  prey  data  and  fine  scale  whale  kinematics  to  be  combined  in  ecological  and  BRS 
analyses. 

•  Initial  statistical  approach  incorporating  echosounder  prey  data  in  models  of  blue  whale 
behavioral  responses  as  an  environmental  covariate. 


RESULTS 

We  conducted  prey  mapping  both  before  and  after  controlled  exposure  experiments  (CEEs)  associated 
with  10  blue  and  2  fin  whales  during  SOCAL-1 1  and  SOCAL-12  that  were  included  in  the  analysis 
depicted  below.  Additionally,  we  have  collected  prey  data  either  before  or  after  behavioral 
response  experiments  on  9  blue  whales. 


Date 

Species 

Pre-exposure 

Post-Exposure 

Control 

7/29/1 1 

Blue  Whale 

X 

X 

7/29/1 1 

Blue  Whale 

X 

7/30/1 1 

Blue  Whale 

X 

X 

7/31/11 

Blue  Whale 

X 

8/1/11 

Blue  Whale 

X 

X 

8/2/11 

Blue  Whale 

X 

8/2/11 

Blue  Whale 

X 

X 

8/3/11 

Blue  Whale 

X 

8/3/11 

Blue  Whale 

X 

8/6/11 

Blue  Whale 

X 

X 

8/6/1 1 

Blue  Whale 

X 

X 

8/7/1 1 

Blue  Whale 

X 

8/8/1 1 

Blue  Whale 

X 

8/9/1 1 

Blue  Whale 

X 

X 

8/9/1 1 

Blue  Whale 

X 
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8/4/12 

Blue  Whale 

X 

X 

8/4/12 

Blue  Whale 

X 

8/4/12 

Fin  Whale 

X 

X 

10/18/12 

Blue  Whale 

X 

X 

10/20/12 

Fin  Whale 

X 

X 

7/26/13 

Blue  Whale 

X 

X 

X 

8/2/13 

Blue  Whale 

X 

X 

X 

8/4/13 

Fin  Whale 

X 

X 

9/14/13 

Fin  Whale 

X 

X 

9/15/13 

Fin  Whale 

X 

X 

X 

9/15/13 

Fin  Whale 

X 

X 

X 

9/16/13 

Fin  Whale 

X 

We  are  presently  in  the  final  stages  of  the  SOCAL-13  field  effort  and  have  collected  additional 
before/after  CEE  prey  measurements  from  2  blue  whales  and  4  fin  whales  as  well  as  one  before  only 
measurement  with  1  fin  whale. 

We  have  explored  the  fine  scale  kinematics  and  energetics  of  blue  whale  feeding  behavior  using  the 
Dtag  sensor  (e.g.  accelerometer)  data  relative  to  acoustically  measured  prey  patches.  Describing  this 
baseline  behavior  is  critical  for  understanding  how  blue  whale  foraging  decisions  vary  as  a  function  of 
prey  distribution  and  abundance.  We  have  found  that  1)  blue  whales  are  targeting  krill  patches  with 
respect  to  patch  density  more  than  patch  depth  building  upon  previous  papers  that  have  shown  a  strong 
depth  selection  independent  of  krill  density  (Croll  et  al.  2002;  Doniol-Valcroze  et  al.  2011).  As  krill 
densities  increase,  the  number  of  feeding  lunges  per  dive  increases  as  well  even  though  other  prey 
patches  are  available  to  a  foraging  whale  (Figure  4).  This  implies  that  understanding  prey  patch  quality 
is  critical  in  assessing  if  whales  conform  to  optimal  foraging  theory  in  the  marine  environment.  In 
addition,  whale  feeding  lunges  were  more  dynamic  in  their  kinematic  motion;  more  likely  to  have 
increased  pitch  and  roll  during  lunges  when  prey  patch  density  was  lower;  higher  densities  of  prey  led 
to  straight  lunges  compared  to  the  acrobat  180°  and  360°  rolls  (whale  and  prey  from  7/29/2011  and 
8/01/2011;  Figure  5). 
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4  3  2  1  100  200  300  400  500 


60  50  40  30  20  10  0  0  5  10  15 

Number  of  Lunges  Ki4  patches  by  depth  {#) 

Figure  4.  a)  Number  of  lunges  and  lunges  per  dive  (blue)  and  b)  krill  patch  depth  and 

krill  patch  density  (red). 


Figure  5.  Dense  patches  of  krill  were  correlated  with  straight  lunges  (above)  compared  to  the  more 
acrobatic  180  and  360°  rolls  at  lower  patch  densities  (below). 
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We  used  generalized  additive  mixed  models  (GAMM)  with  individual  whale  as  a  random  covariate  to 
assess  the  effect  of  prey  density  and  patch  structure  on  whale  feeding  kinematics.  Results  from  pitch 
and  roll  models  indicate  that  whales  are  more  acrobatic  when  foraging  on  less  dense  patches  (Figure  6 
below).  All  of  this  is  critical  in  putting  behavioral  responses  in  context  to  the  baseline  ecology  of  these 
animals. 

We  have  completed  analyses  on  the  prey  mediated  response  to  sonar  playbacks  using  the  PCA-GAMM 
approaches  described  in  Goldbogen  et  al.  2013.  Given  we  have  prey  data  only  before  and  after  the 
playbacks,  we  were  not  able  to  include  a  “during”  state  so  any  behavioral  state  changes  that  returned  to 
normal  upon  cessation  of  playback  would  not  be  detected.  The  GAMM  examining  response  in  dive 
axis  1  (dive  time,  surface  time,  breaths,  dive  depth,  etc.)  showed  a  significant  before-after  effect 
including  potential  changes  in  prey.  However,  we  found  that  bottom  depth,  school  height,  and  school 
depth  explained  significantly  more  variability  in  dive  axis  1  than  playback  alone  (r2  =  0.703  for  the 
best-fit  model  shown  in  Table  1,  Figure  7);  treatment  type  (MFA  or  PRN)  was  not  significant  for  any 
of  the  models.  This  makes  the  argument  that  a)  prey  are  critical  in  understanding  behavioral  changes 
(r2  with  prey  of  0.703  compared  to  r2  of  0.14  for  PCA-GAMM  on  Dive  Axis  1  in  Goldbogen  et  al. 
2013)  and  b)  we  prey  measurements  during  exposure  would  be  useful  to  collect  as  many  whales 
returned  to  pre-exposure  behavior  immediately  following  acoustic  playbacks.  Given  the  high  amount 
of  variance  explained  by  the  mathematical  models  that  include  prey  parameters,  we  will  also  examine 
how  prey  changes  may  have  affected  individual  whale  responses  in  the  Mahalanobis  distance  /  change- 
point  analysis  approach  currently  underway  by  the  SoCal  BRS  team.  As  these  techniques  are  still  being 
developed,  we  have  prepared  the  prey  framework  for  inclusion  but  are  waiting  for  the  complete 
analyses.  Together  the  PCA-GAMM  across  individual  analyses  and  the  individual  time  series  analyses 
will  comprise  a  complete  manuscript  targeted  to  a  journal  such  as  Proceedings  B  Biology  for 
submission  in  late  2013/early  2014. 


Nurtoer  of  Patches  (tf/10  rrwi) 


Figure  6.  Blue  Whale  Response  mediated  by  prey 
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Table  1.  GAMM  summary  for  best-fit  model  with  Dive  Axis  1  as  a  function  of prey 

and  playback  state. 


Degrees  of 
freedom 

F  statistic 

p-value 

Before-After 

1.1803 

3.271 

0.0017 

** 

s(Sv mean) 

4.95E-07 

0 

0.388375 

s(Height mean) 

9.40E-01 

3.501 

0.000147 

*** 

s(Depthmean) 

2.47E+00 

17.595 

1.68E-12 

*** 

s(BotDep) 

9.79E-01 

9.848 

8.91E-09 

*** 

r2=  0.703 

Svjnean 


Height_mean 


Figure  7.  Dive  axis  1  generalized  additive  mixed  model  output  showing  how  dive  behavior  changes 

significantly  as  prey  patches  change. 


IMPACT/IMPLICATIONS 

Our  research  to  date  has  significantly  increased  our  understanding  of  baleen  whale  feeding  behavior 
and  foraging  ecology  (Friedlaender  et  al.  in  review,  Hazen  et  al.  in  prep). 

We  have  also  conducted  initial  analyses  that  will  contribute  significantly  to  our  ability  to  understand 
how  the  fine-scale  kinematics  of  baleen  whale  feeding  are  affected  by  changes  in  the  distribution, 
abundance,  and  behavior  of  their  prey.  This  information  will  provide  critical  baseline  knowledge  to 
interpret  behavioral  responses  of  baleen  whales  to  anthropogenic  sounds.  Understanding  how  baleen 
whale  foraging  conforms  to  ecological  theory  (optimal  foraging  theory)  and  how  changes  in  prey 
behavior  affect  the  complexity  of  feeding  behaviors  will  increase  our  understanding  of  baleen 
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whales  significantly  can  be  compared  to  similar  field  studies  across  other  species  and  in  different 
locations. 

Finally,  our  latest  analyses  augment  our  initial  findings  of  the  behavioral  responses  of  blue  whales  to 
anthropogenic  sounds.  Specifically,  including  quantitative  metrics  of  prey  patches  being  targeted  by 
whales  actively  feeding  during  behavioral  response  experiments  significantly  increases  the  predictive 
power  of  our  statistical  models.  In  our  recently  published  blue  whale  response  paper  (Goldbogen  et  al. 
2013),  the  behavioral  response  analysis  including  whale  kinematic  data  but  no  prey  data  had  an  r2  of 
0.14.  Running  the  same  model  (GAMM)  including  prey  patch  metrics  increased  the  r2  to  0.703,  and 
the  most  significant  factors  in  this  model  were  all  related  to  prey.  Response  to  playback  was  still 
significant  given  prey  covariates,  which  corroborates  Goldbogen  et  al.  (2013),  but  specific  examination 
of  individual  responses  is  critical  to  further  tease  this  apart.  Thus,  we  have  achieved  greater 
understanding  of  prey-mediated  behavioral  responses  in  baleen  whales.  This  is  important  information 
for  the  behavioral  response  study  because  with  such  increased  explanatory  power  for  when  changes  in 
whale  behavior  can  be  detected,  we  can  more  specifically  characterize  how  whale  behavior  is  affected 
by  sound  exposure. 

Our  results  relating  to  how  the  fine-scale  kinematics  of  baleen  whale  feeding  are  affected  by  prey 
behavior  will  be  included  in  a  manuscript  to  be  submitted  in  a  high  profile  journal  such  as  the 
Proceedings  of  the  National  Academy  of  Sciences  in  late  2013. 

RELATED  PROJECTS 

This  project  is  closely  coordinated  with  the  Southern  California  Behavioral  Response  Study  (SOCAL- 
BRS  -  see:  www.socal-brs.org)  which  is  measuring  behavior  and  responses  to  simulated  mid¬ 
frequency  sonar  and  other  signals  in  marine  mammals.  The  prey-mapping  measurements  for  this 
project  are  leveraging  boat  time  and  other  logistical  support  from  the  ongoing  SOCAL-BRS  project, 
while  providing  data  that  are  directly  relevant  to  interpreting  the  behavioral  responses  of  mysticetes  to 
controlled  exposure  experiments. 
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